Nonlinear optical response in two-dimensional Mott insulators 



Makoto Takahashi, Takami Tohyama and Sadamichi Maekawa 
Institute for Materials Reserch, Tohoku University, Sendai 980-8577, Japan 
(Dated: February 6, 2008) 

We study the third-order nonhnear optical susceptibihty x'^' ^'iid photoexcited states of two- 
dimensional (2D) Mott insulators by using an effective model in the strong-coupling limit of a 
half-filled Hubbard model. In the numerically exact diagonalization calculations on finite-size clus- 
ters, we find that the coupling of charge and spin degrees of freedom plays a crucial role in the 
distribution of the dipole-allowed states with odd parity and the dipole-forbidden states with even 
parity in the photoexcited states. This is in contrast with the photoexcited states in one dimension, 
where the charge and spin degrees of freedom are decoupled. In the third-harmonic generation 
(THG) spectrum, the main contribution is found to come from the process of three-photon reso- 
nance associated with the odd-parity states. As a result, the two-photon resonance process is less 
pronounced in the THG spectrum. The calculated THG spectrum is compared with recent experi- 
mental data. We also find that x'^"* with a cross-polarized configuration of pump and probe photons 
shows spectral distributions similar to x'"^' with a copolarized configuration, although the weight is 
small. These findings will help analyses of the experimental data of x'^' in 2D Mott insulators. 
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I. INTRODUCTION 

The charge gap in Mott insulators is a consequence of a 
strong electron correlation represented by a large on-site 
Coulomb interaction. The correlation induces novel phe- 
nomena in terms of the interplay of charge and spin de- 
grees of freedom. 0] In one-dimensional (ID) Mott insu- 
lators, two particles created by photoexcitation — i.e., an 
unoccupied site and a doubly occupied site of electrons — 
can move inside the system without being disturbed by 
surrounding spins in the background. This is a mani- 
festation of a separation of the charge and spin degrees 
of freedom, called the spin-charge separation inherent 
in ID correlated electron systems. Optical responses in 
the ID Mott insulators are characterized by this phe- 
nomenon. ||, |, |, |, |, In two-dimensional (2D) Mott 
insulators, on the other hand, the two particles are ex- 
pected not to be free from the spin degree of freedom, 
because the propagation of a carrier is known to induce 
a spin cloud around the carrier as a consequence of the 
misaligned spins along the carrier-hopping paths. Q Such 
an interplay of spin and charge is also one of main sub- 
jects of the study of high-temperature superconductivity. 

The nature of the two particles in the photoexcited 
states of the Mott insulators is obtained by examining 
the linear susceptibility x''^^ with respect to the applied 
electric field, which provides information on the dipole- 
allowed states with odd parity among the photoexcited 
states. In addition to x^^\ the third-order nonlinear opti- 
cal susceptibility x^^^ is useful to detect not only the odd- 
parity states but also the dipole-forbidden states with 
even parity. |^ Recently, large values of x*-'^^ have been 
reported for ID Mott insulators of copper oxides such as 
Sr2Cu03 ( Ref. 9 and 10 ) from an electro-reflectance 
measurement 0] and pump and probe spectroscopy. 
Analyses of x^^bave suggested that odd- and even-parity 
states are nearly degenerate with a large transition dipole 



moment between them. Theoretically, x'-"^-' in ID Mott 
insulators has been examined by employing the numeri- 
cally exact diagonalization technique for small clusters of 
the Hubbard model at half filling. |6| It has been shown 
that odd- and even-parity states are almost degenerate 
in the same energy region and that the degeneracy is due 
to the spin-charge separation and strong on-site Coulomb 
interaction. [pd]| From this viewpoint, an effective model 
that can describe the optical nonlinearity has been pro- 
posed, |6| that is, a holon-doublon model, where holon 
and doublon represent the charge degree of freedom for 
photoinduced unoccupied and doubly occupied sites, re- 
spectively. The model reproduces very well the charac- 
teristic behaviors of the experimental x''^'' including data 
from third-harmonic generation (THG) spectroscopy. 0| 

In the 2D Mott insulators of copper oxide such as 
Sr2Cu02Cl2, x^^^ has been reported to be one order of 
magnitude smaller than that in ID. ||l^ Theoretical cal- 
culations based on the (extended) Hubbard model have 
shown such a dimensionality dependence of x''^'' ■ 
However, a complete understanding of the nonlinear op- 
tical responses in 2D Mott insulators has not been ob- 
tained both experimentally and theoretically. For in- 
stance, THG spectra in 2D have been reported from three 
groups, |l^, |lj, |l^ but the data have not converged: 
One of the data in La2Cu04 shows several pronounced 
peak structures in the energy region of 0.6 - 1.1 eV,|l^ 
while other data in La2Cu04 (Ref. 14) and Sr2Cu02Cl2 
(Ref. 15) exhibit a broad structure with a maximum at 
around 0.7 eV. Although the former data of THG were 
analyzed by using an excitonic cluster model |^6j taking 
into account both copper and oxygen ions, a microscopic 
description and understanding of photoexcited states in 
2D are not complete in terms of the interplay of charge 
and spin degrees of freedom. 

In this paper, we theoretically examine photoexcited 
states and nonlinear optical responses in the 2D Mott 
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insulators. We use an effective Hamiltonian of a half- 
filled Hubbard model in the strong-coupling limit and a 
numerically exact diagonalization method on finite-size 
clusters for the calculation of x^^"^ ■ By using the effective 
Hamiltnian, we can treat the clusters larger in size than 
those used in a previous study for the 2D Hubbard model 
in Ref. 4, leading to new insights into the photoexcited 
states in the 2D Mott insulators. It is found that the 
edge of the distribution of the even-parity states is lo- 
cated lower in energy than that of the odd-parity states. 
This is different from the distribution in ID, where the 
edges of both states almost coincide. The origin is at- 
tributed to the presence of an exchange interaction be- 
tween spins. In other words, the spin degree of freedom 
plays an important role in the photoexcited states. In 
the THG spectrum, dominant contributions come from 
the process of three-photon resonance associated with the 
odd-parity states. The two-photon resonance process is 
hidden by the dominant contributions. The spectrum ob- 
tained by using realistic parameters for 2D copper oxides 
thus shows broad maxima coming from three-photon res- 
onance. Such broad spectral features are qualitatively in 
agreement with the experimental data in Refs. 12 and 
13. We also demonstrate that x^^^ obtained by a cross- 
polarized configuration of pump and probe beams shows 
spectral distributions similar to x*-"^) with a copolarized 
configuration but with small weight. These findings will 
help the analysis of the experimental data of x^'^) in the 
2D Mott insulators. 

The rest of this paper is organized as follows. We intro- 
duce an effective Hamiltonian of the half-filled Hubbard 
model in the strong-coupling limit and show outlines of 
the procedure to calculate x^^^ X^^'' in Sec. II. In 
Sec. HI, calculated results of the linear and nonlinear op- 
tical responses are presented. The distributions of pho- 
toexcited states are discussed in terms of the effect of the 
spin degree of freedom. Two-photon absorption (TPA) 
and THG spectra are calculated with both copolarized 
and cross-polarized configurations. We also compare our 
results with existing experimental data of x^^^ • The sum- 
mary is given in Sec. IV. 

II. MODEL AND METHOD 

The insulating cuprates are known to be charge- 
transfer- (CT-) type Mott insulators, where both 3d and 
2p orbitals exist in the transition- metal and ligand ions, 
respectively, and participate in the electronic states. The 
value of the gap is predominantly determined by the en- 
ergy position of the 2p orbitals. However, it is well es- 
tablished that the electronic states of the CT-type insu- 
lators can be described by an extended Hubbard model 
with a half-filled single band by mapping a bound state 
called the Zhang-Rice singlet state onto the lower Hub- 



bard band. Q The extended Hubbard Hamiltonian in 2D 
is given by 

E {4,aCj,a + ii-C.)+V J2 n,nj, (1) 

where ^ is the creation operator of an electron with spin 

a at site i, n^^a = c\^^Ci^a and ni=ni^^+ni^i, {i,j)ist runs 
over pairs on the nearest-neighbor (NN) sites, 2nd 
runs over pairs on the next-nearest-neighbor (NNN) sites, 
t {t') is the NN (NNN) hopping integral, U is the on-site 
Coulomb interaction, and V is the Coulomb interaction 
between the NN sites. The value of t is estimated to be 
t ~ 0.35 eV from analysis of the electronic structures in 
the cuprates. t' has been systematically found to be 
of negative sign in contrast to the NN hopping ampli- 
tude t with positive sign and to have 40% of its magni- 
tude {t'/t = -0.4). H The value of U is estimated to be 
U —lot for the gap values to be consistent with experi- 
mental ones. V/t is taken to be 1. [0 

In the strong-coupling limit {U ^ t), the ground state 
at half filling has one spin per site; i.e., there is no doubly 
occupied site. In this case, the low-energy excitation of 
a Hubbard model may be described by the Heisenberg 
model, where the dominant exchange interaction is the 
NN interaction J given by At'^/U. The Hamiltonian reads 




where Si is the spin operator with S = 1/2 at site i. 
However, the photoexcited states that are created by the 
light have both one doubly occupied site and one va- 
cant site. In order to obtain an effective Hamiltonian of 
the extended Hubbard model that describes the photoex- 
cited states, we introduce projection operators Hq, Hi, 
and II2 onto the Hilbert space with no doubly occupied 
site, one doubly occupied site, and two doubly occupied 
sites, respectively. The Hubbard Hamiltonian, Eq. (|l|), is 
conveniently split into two parts: kinetic terms (Ht and 
Ht') and interacting terms {Hu and Hy)- We perform 
a second-order perturbation with respect to Ht- Since 
the value of \t'\ is smaller than t {t' = —OAt), we neglect 
Ht' in the perturbation process. The resulting effective 
Hamiltonian is given by 

i/eff = Hi+H2 + H3 + Hi + U, (3) 

where 



3 



Ili{Ht+Ht,)Ili 



1 

c7 



(i,i)2nd:0- 



'jo-} 



i.j,k,a.(T' 

+ (1 - nj_„)c'-Ci^ni_„nk-a'cl^,Cj„,{l - Uj^^,)}, 



t 



Ckcnk- 



(4) 



(5) 



+nicrc|_^4^(l - nk^a)CjaCj-a - (1 - nk-a)c}f.^Ckank-a{l - ni-a)nj-anja}, 



(6) 



and 



iJ4 = -V"^n,|n,x(l -njt)(l -rijx). 



(7) 



Here, 77i expresses hopping process of unoccupied and 
doubly occupied sites. H2 contains the exchange inter- 
action and so-called three-site terms. Hi and H2 are 
similar to the t-J model with three-site terms and NNN 
hopping, while Ht, has specific terms of the photoexcited 
states. The attractive interaction between the unoccu- 
pied and doubly occupied sites is described by H4. 

The electric field E of the incident light induces the di- 
electric polarization P in a material, which is described 
by a power series of nonlinear optical susceptibility x*^") : 
P = eo(x^^^E + x^^^E^ + x'-^'^E^ H ). The linear sus- 
ceptibility x^^^ is given by 
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where eg is the dielectric constant, k and / are the po- 
larization directions and eroa{— e(0|r|a)) is the dipole 
moment between the ground state |0) of the Heisenberg 
model (^) and photoexcited state |a) with odd parity ob- 
tained from the effective Hamiltonian (^), f being the 
dipole displacement operator. fla{^ Ea — Eq) is the 
energy difference between |0) and |a), L is the number 
of sites, and Fa is the damping factor. The distribu- 
tion of odd-parity states among the photoexcited states 
is obtained by examining the imaginary part of x^u'^ in 



the / polarization direction, which is related to the dy- 
namical current-current correlation function au{oj) as Im 
Xu\-^]^) = oiii{u!) I u? . Here an(w) reads 



where j; is the / component of the current operator j 
which is given by 

j = E (I^J "I^»)(4<TCi,<T 

Here, Ri is the position vector at site i. The distribu- 
tion of even-parity states can be detected by a correlation 
function fiuiiS) in which the current operator jj in auiiS) 
is replaced by the / component of the stress tensor oper- 
ator f , 



/3nH = E K^l^'|0)l''5(^ - ^6 + ^0), (11) 

h 



with 
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where 6 denotes even-parity state |6) with energy . Due insulating cuprates belong. The third-order susceptibil- 
to symmetry restrictions, the second-order susceptibility ity x'^'^^ is expressed as 
X^^) vanishes in centrosymmetric materials to which the 
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where c denotes odd-parity states |c) with energy Ec, 
uja = uJi + 0J2 + ^z, and V represents the sum of per- 
mutation on {l,uji), (171,102), (n^uj^), and (k,u!a-)- Here- 
after, the damping factors Fq, Tb and Tc are assumed 
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to have the same value F, with T/t — 0.4, and we take 
e = ft = eo = 1. 
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in Eq. (nih is rewritten as 
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(15) 



where Zp = u}p + iV. This quantity can be calculated by using the correction vector technique [|l8[ : 

(0| f \a) {a\ f' |5) (5| f" |c) (c| f" |0) 
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where r'^ is the fc component of f . The correction vectors 
(j){zr) and i^izq, Zr) are obtained by solving the following 
equations: 

{H,f,-Ea~Zr)\<j){Zr)) = P|0) (17) 
{Hcfi-Eo~Zq)\lP{Zq,Zr)) = f"|(/.(z,)), (18) 

where f"|0) is obtained by solving j„|0) = i(i?cff — 
So)f"|0). 

In order to describe 2D systems, we employ an 18-site 
cluster with open boundary condition, which is shown 
in the inset of Fig. 1(a). The ground state |0) of this 
cluster is obtained by applying the standard Lanczos 
method to the Heisenberg Hamiltonian (||). The cor- 
relation functions aii{Lij) and Pii{ui) are also calculated 
by the Lanczos technique, x^'^^ is obtained from Eqs. 
(p^)-(p^. In the following section, we consider two 



kinds of x^'^^ spectra: One is theTPA spectrum which 
is defined as Imx'^' (— w; — w, w), and the other is 
THG spectrum which is defined as |x^^H~3a;; w, w, w)|. 
Since the 18-site cluster has no 90° rotational symmetry, 
two spectra in the x- and y-polarization directions are 
inequivalent. However, qualitatevely similar features 
are obtained between the two spectra. We have also 
calculated the x^'^^ spectra by using a 4x4 cluster that 
preserves the 90° rotational symmetry. The results 
are found to be consistent with those of the 18-site 
cluster, indicating small size effects. Therefore, in 
the next section, we take an average of the spectra 
of the 18-site cluster and use the following notation: 

Inix(i) (-w; uj) ={ImxS(-tJ; uj)+lmx^y^ t^)}/2, 
a{uj) = {a^xii^) + ayy{uo)]/2, /3{uj) = {/3xx{(^) + 
/5aa(w)}/2, xj^^ = {lmx''xxxx{-i^; -io,uj,uj)+ 
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+ IXTOU(-3w;a;,a;,w)|}/2. 
cross-polarized configurations of the pump 
probe photons, we take the following av- 
erages: Xl^"^ = {lmX^xxyyi~(^;~UJ,UJ,Uj) + 

^T^Xyyxx{-u};-uj,uj,uj)}/2 and = 
{|xS'yy(-3w;w,w,w)| -|- |xyyL(-3w; w, w, w)l}/2. 



III. RESULTS AND DISCUSSIONS 

A. Linear absorption and two-photon absorption 

We first show the distribution of optical-allowed states 
with odd parity and optical-forbidden states with even 
parity among the photocxcitcd states created by pho- 
toexcitation. Figure 1(a) exhibits a{o-!), which detects 
the odd-parity states. Here 6 functions denoted by ver- 
tical bars are broadened by a Lorentzian with a width 
of OAt. There are two broad peaks a.t lu ~ U ^ —3t 
and —t, which is in contrast with ID where only a single 
peak appears [see Fig. 3(b)]. We note that the higher- 
energy peak plays important roles in the resonance condi- 
tion of two-magnon Raman scattering as will be reported 
elsewhere. [|l9| The distribution of even-parity states de- 
tected by (3{uj) is shown in Fig. 1(b). The lowest-energy 
broad peak is located at w — C/ ~ —3.5t, which is lower 
than that of a{uj). This is also different from the case of 
ID where a one-to-one correspondence between the odd- 
and even-parity eigenstates is clearly seen. [|[ 0, |ll| Such 
a degeneracy in ID is a consequence of the spin-charge 
separation and strong on-site Coulomb interaction. In 
2D, however, motion of charge carriers induces a spin 
cloud around the carriers as a consequence of the mis- 
aligned spins along the carrier-hopping paths. Therefore, 
the effect of the coupling of charge and spin should ap- 
pear in the distribution of odd- and even-parity states. 
This effect can be checked by changing the value of the 
exchange interaction J/t. In Fig. 2, we show a{Lu) and 
(3(u}) with the same parameters as used in Fig. 1 but 
U/t = 100 {J/t = 0.04). By changing J/t from 0.4 to 
0.04, we find that the higher-energy peak at uj — U = — t 
in Fig. 1(a) is smoothly connected to a broad peak at 
UJ — U — ~2.5t in Fig. 2(a). On the other hand, the 
lowest-energy peak in Fig. 1(a) loses its weight with re- 
ducing J/t, and only a small amount of the weight is 
seen around lo — U = —bt in Fig. 2(a). Similarly, the 
lowest-energy peak of (3{lo) seen in Fig. 1(b) almost dis- 
appears in Fig. 2(b). Therefore, for small J/t, it is hard 
to identify the difference in the distribution of the odd- 
and even-parity states near the edge of the photoexcited 
states. This implies that the distribution near the edge 
is sensitive to the coupling of charge and spin. In other 
words, the coupling could be the origin of the lowest- 
energy peaks in Fig. 1 and also of the lower-energy shift 
of the even-parity states as compared with the odd-parity 
states. 
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FIG. 1: (a) aiuj) and (b) I3{ijj) in an 18-site cluster of an 
effective model in the strong-coupling limit of the Hubbard 
model with U/t = 10 {J/t = 0.4), t' /t = -0.4, and V/t = 
1. The solid lines are obtained by performing a Lorentzian 
broadening with a width of QAt on the 5 functions denoted 
by vertical bars. The spectra in (a) and (b) represent the 
distribution of photoexcited states with odd and even parities, 
respectively. The inset in (a) shows the 18-site cluster with 
open boundary condition. 



The pump and probe spectroscopy is employed to de- 
tect even-parity states. |l^, ^ As a corresponding quan- 
tity, we calculate the TPA spectrum xj[^^- In ^\^^^ ^ 
dominant contribution comes from the following term: 
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(f^a — iF — ijj){U,b — iF — 2w)(ric — iF ■ 



(19) 

where the dipole moment between |0) and \h) in the x{y) 
direction is represented by XQa{yoa) instead of ?'oa(''oa) 
in Eq. (^3|). The resonance occurs when 2t<j is equal to 
the energy of eigenstates with even parity. Figure 3(a) 
shows xj^^ (solid line) together with linear absorption 

Im x^^-* (~2a;; 2a;) (dashed line). In the figure, the two 
spectra are plotted by making the largest weight the same 
height. In order to examine the effect of dimensionality, 
results for an 18-site ID chain are also shown in Fig. 3(b), 
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FIG. 2: The same as Fig. 1, but U/t = 100 {J/t = 0.04). 



where parameter values are those of the ID Mott insu- 
lator SraCuOg (C//t=10, V lt=l.b, and t' /t = 0). The 
maximum value of TPA in 2D is nearly 5 times smaller 
than that in ID. We also find that, while in ID a peak 
position as well as an edge position is nearly identical 
between linear absorption and TPA spectra, in 2D both 
peak and edge positions are lower in TPA than in linear 
absorption as expected from Fig. 1. Such dimensionality 
dependence is qualitatively consistent with experimental 
data where a low-energy tail of TPA in Sr2Cu02Cl2 down 
to 1.5 eV was reported. Q 

Summarizing the dimensionality dependence of optical 
responses in the Mott insulators, we find that (i) the lin- 
ear absorption spectrum in 2D shows two broad peaks 
in contrast to a single peak in ID, (ii) the magnitude of 
the TPA spectrum in 2D is nearly 5 times smaller than 
that in ID, and (iii) in 2D the peak and edge positions 
of the TPA spectrum are lower than those of linear ab- 
sorption in contrast to the case of ID, where TPA and 
linear absorption show almost an identical behavior near 
the edge. All of the features are consistent with the ex- 
perimental data of linear absorption and TPA spectra 
measured in ID SraCuOg and 2D SrzCuOzCla. |l2| This 
agreement implies that our single-band Hubbard model 
is suitable for a description of the photoexcited states in 
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FIG. 3; (a) xf^^ (solid line), xF^ (dotted line), and 

Imx*^'(-2t^;2a;) (dashed line) in an 18-site 2D cluster of an 
effective model in the strong-coupling limit of the Hubbard 
model with V jt = 10, t'/t = -0.4, and Vjt = 1. The 
solid and dashed arrows denote the peak position of Xy"^* 
and Imx^^"* (— 2tj; 2a;), respectively, (b) Xy"^^ (solid line) and 
Imx'^H-2Lj;2a;) (dashed line) in an 18-site ID chain with 
Ult = 10, t' It = 0, and Vjt = 1.5. 



both ID and 2D Mott insulators of copper oxides. In the 
next subsection, we will examine THG spectra in 2D, for 
which three experimental results are now available but 
exhibit some discrepancies between them. 



B. Third-harmonic generation 

The THG spectrum x™*^ contains information about 
both odd- and even-parity states via multiphoton res- 
onance processes, since a dominant contribution is ex- 
pressed as 

1 Xi)aXabXhcXc<^ yQayabybcUcO 

2L {^a -iV- 3tj)(f76 -iV- 2uj){Vlc -iV-ujy 

(20) 

where the resonance occurs when Sw (three-photon res- 
onance) and 2uj (two-photon resonance) are equal to the 
eigenenergies of the odd- and even-parity states, respec- 
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FIG. 4: x™'^ (solid line) and x™'^ (dashed line) in an 18-site 
2D cluster of an effective model in the strong-coupling limit of 
the Hubbard model with U/t = 10, t' /t = -0.4, and V/t = 1. 
The dashed arrows denote the position of three-photon reso- 
nance expected from Imx^^'' {—'2uj;2lu) in Fig. 3(a). The solid 
arrow denotes the position of two-photon resonance expected 
form xj^^- 



tively. Figure 4 shows x™^ (solid line). The spectrum 
consists of a broad peak and a hump below u ^ 3.5t. 
Since the peak position of a; 2.3t (dashed arrow) coin- 
cides with one-third of a peak energy at w ~ 6.9t in the 
hnear absorption spectrum indicated by a dashed arrow 
in Fig. 3(a), the origin of the peak can be attributed to 
three-photon resonance. On the other hand, there are 
two possibilities for the origin of the broad hump struc- 
ture around lu ~ 2.8t. One is again three-photon reso- 
nance, and the other is two-photon resonance resonating 
with even-parity states. The dashed arrow at u — 2.9t 
in Fig. 4 indicates one-third of the second-peak energy 
in the hnear absorption spectrum (8.8t) of Fig. 3(a) indi- 
cated by a dashed arrow, while the solid arrow at a; = 3.2t 
indicates the energy of the lowest-energy peak in the 
TPA spectrum of Fig. 3(a) indicated by a solid arrow. 
Since the hump position is closer to the dashed arrow, 
we regard the hump as a result of three-photon reso- 
nance. A contribution of the even-parity states through 
two-photon resonance is supposed to overlap with three- 
photon resonance and thus may be hidden as a back- 
ground in the THG spectrum. 

We can check the contribution from the even-parity 
states to THG, by examining the real part of x™*^: 
which is shown in Fig. 5(a). The solid curve obtained 
from Eq. ([l^ ) changes its sign from negative to positive 
at a; ~ 3.2i, which is close to the two-photon resonance 
energy (the solid arrow in Fig. 4). The contribution from 
the dominant term, Eq. ( pO| ) is also plotted as the dashed 
curve, which is similar to the solid curve and shows a zero 
crossing at lu ^ 3.3i . In order to clarify whether such 
a crossing or sign change is due to a contribution from 
the even-parity states, we calculate the real part of THG 
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FIG. 5: The real part (a) and the phase (b) of X™'' in ^-n 
18-site 2D cluster of an effective model in the strong-coupling 
limit of the Hubbard model with U/t = 10, t' /t = -0.4, and 
V/t — 1. The solid lines represent full contribution obtained 
from Eq. (p^. The dashed lines represent the contribution 
from the dominant term, Eq. ([2o|). The dotted lines are 
obtained by neglectinethe two-photon resonance terms in the 
denominator of Eq. (EOl). 



by neglecting the two-photon resonance contribution in 
the denominator of Eq. (|2^) and plot it as the dotted 
curve in Fig. 5(a). We find no sign change near the two- 
photon resonance energy. Therefore, the sign change is a 
consequence of the presence of the even-parity states in 
the THG spectrum. However, as mentioned above, the 
weight cannot be recognized as a pronounced structure 
because of large contributions from the odd-parity states 
in the same energy region. 

The contribution from the even-parity states to the 
THG spectrum can also be seen through the phase 9 of 
the THG susceptibility: x^^^ = |x^^^|e'^- In Fig. 5(b), the 
solid line represents obtained from the data of Eq. (|l^) . 
Here 9 has two hump structures at w ~ 2.3i and 3.2i. The 
hump structure at (jj ~ 2.3i is caused by a three-photon 
resonance associated with the peak at the same energy 
in x™^- The position of the other hump structure at 
LU ^ 3.2t is close to a two- photon resonance energy. The 
phase obtained from the dominant term, Eq. ([20|), is also 
plotted by the broken curve, which is similar to the solid 
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curve. The higher-energy hump structure disappears (as 
the dotted curve) the contribution from two-photon res- 
onance is neglected in the denominator of Eq. (EG) . 



C. Cross-polarized TPA and THG 

The 18-site cluster employed here has two mirror 
planes associated with the x and y axes [see the inset 
of Fig. 1(a)] and thus have a C2v point group whose ir- 
reducible representation is Ai, A2, Bi, and -B2. Since 
the ground state of the Heisenberg model |0) belongs to 
Ai, the photoexcited state with odd parity \a) appearing 
in the dipole moment xoa ivoa) belongs to Bi {B2). On 
the other hand, the states with even parity \b) encoun- 
tered in x™^ and xj[^^ belong to Ai. In other words, 
the eigenstates of Bi and B2 representations contribute 
to both one- and three-photon resonance, and those of 
Ai emerge as two-photon resonance. Therefore, there 
is no information on A2 subspace in xj^'^ ^^'^ X™*^- 
However, if we investigate cross-polarized nonlinear sus- 
ceptibilities Xa^^ ^iid x™*^; are able to extract in- 
formation about A2 eigenstates because the state \b) be- 
longs to either Ai or ^2- In xl^^: the numerator of 
the dominant term, Eq. (|l9|), becomes XQayab{xbcycO + 
ytcXco) + yoaXabiybcXcO + Xbcyco), leading to two-photon 
resonance associated with A2 subspace at ilf, = 2lu. In 
contrast, the numerator of the dominant term of x™*^ 
[Eq. (|^)] is the sum of xoaXabybcycO + yoayabXbcXcO and 

XQayabixbcyco + ybcXco) + yoaXabiybcXco + Xbcyco), where 

two-photon resonance occurs within Ai and A2 sub- 
spaces, respectively. 

In Fig. 3(a), xT^'^ is plotted as the dotted curve. As 
mentioned above, almost all of weight comes from the 
A2 representation. The weight is smaller than that of 



TPA 



and a maximum appears in the higher-energy re- 
gion around to = At, although the weight spreads over in 
the energy range similar to X|[^^- The x™*^ spectrum 
is shown in Fig. 4 as the dashed curve. Similar to X™'~^i 
the contribution from three-photon resonance is domi- 
nated as evidenced by the energy position of the weight. 
The Ai contribution of two-photon resonance overlaps 
the three-photon resonance contribution and thus is hid- 
den as a background. The A2 contribution, on the other 
hand, should appear around lu = At according to the dis- 
tribution in TPA, but the weight is not enhanced. The 
total weight is smaller than that of x™^ i being similar 
to the TPA spectra. In order to obtain further insight 
into the photoexcited states of the 2D Mott insulators, 
it is desired to measure experimentally the weight and 
distribution of the cross-polarized TPA and THG and 
compare them with our results. 

IV. SUMMARY 

We have examined the photoexcited states and nonlin- 
ear optical responses in the 2D Mott insulators by using 



an effective model for a single-band Hubbard model in the 
strong-coupling limit. By employing the numerical tech- 
nique for small clusters, we have obtained characteristic 
features that are consistent with the experimental data 
of linear absorption and TPA for 2D insulating cuprates. 
They are the following: (i) linear absorption in 2D shows 
two broad peaks in contrast to a single peak in ID, (ii) 
the magnitude of the TPA spectrum in 2D is smaller than 
that in ID, and (iii) in 2D the peak and edge positions of 
the TPA spectrum are lower than those of linear absorp- 
tion in contrast to ID where TPA and linear absorption 
show almost an identical behavior near the edge. The 
interplay of charge and spin degrees of freedom plays an 
essential role in such dimensionality dependence: In the 
photoexcited states of the ID Mott insulator there is a 
separation of the two degrees of freedom, while in 2D the 
charge motion is strongly influenced by the presence of 
spins in the background. Being consistent with this pic- 
ture, we found that the features (i) and (iii) are sensitive 
to the value of the exchange interactions between spins. 

In the THG spectrum of the 2D Mott insulators, dom- 
inant contributions are found to come from the process 
of three-photon resonance associated with the odd-parity 
states. The two-photon resonance process is thus hidden 
by the dominant contributions. The spectrum obtained 
by using realistic parameters for 2D insulating cuprates 
thus shows broad spectral features, which are qualita- 
tively in agreement with the experimental data with a 
broad maximum, |lj, |l5) but not to the data with pro- 
nounced peak structures. [ [T^ We hope further the exper- 
imental efforts to clarify the characteristic of the THG 
spectrum in 2D. 

We have calculated x'^' under the cross-polarized con- 
figuration of pump and probe beams. It is shown in both 
the TPA and THG spectra that the spectral distribu- 
tion is similar to x'"'^-' with copolarized configuration but 
the weight is small. Measurements of x*-'^'' with the cross- 
polarized configurations are highly desired to obtain deep 
insight into the nature of photoexcited states in 2D Mott 
insulators. 
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